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Abstract 

An allozyme comparison of eight taxa in the genus Gazella (Bovidae: Artiodactyla) was conducted to 
clarify the systematic relationships of endangered gazelles currently bred in Saudi Arabia for 
reintroduction. Electrophoretic variation at 16 genetic loci suggested that several similar taxa of 
Arabian gazelles, namely G. gazella gazella, G. gazella erlangeri, G. gazella farasani , and G. gazella 
cora, belong to the same species of G. gazella sensu lato. Four other species proved to have diverged 
genetically: G. thomsoni and G. dorcas , which cluster together, G. rufifrons, and G. subgutturosa. The 
subgenus Trachelocele , in which the latter species has been placed according to morphological 
characteristics, is not supported. Polymorphism and heterozygosity values found in Gazella were 
generally similar to average values reported for mammals. The results are discussed in terms of the 
strategy to follow a conservation program that take genetic data into account. 


Introduction 

There is urgent need for conservation action concerning gazelles (Ryder 1987), both in 
captive and wild populations. Of the 12 species of Gazella (Corbet and Hill 1980), nine 
are considered vulnerable or endangered (IUCN 1988), mainly due to overhunting and 
habitat destruction. In Saudi Arabia, where at least three "good biological species" (Mayr 
1963) are found: G. saudiya , G. subgutturosa and G. gazella (see Fig. 1), the situation is of 
particular concern (Thouless et al. 1991). Over the last few years, tremendous efforts have 
been undertaken towards the conservation of gazelles in this country (Abu-Zinada et al. 
1989). It is widely accepted that species conservation must be based on proper systematics 
of the endangered taxa (Ryder 1986; Groves 1988). 

Morphological characters such as size and shape of horns and skull have been previ¬ 
ously used to establish systematics in Gazella (Groves 1969, 1983; Lange 1972; Groves 
and Lay 1985; Alados 1986/1987). Color and coat patterns have also been used, but these 
characters may vary due to environmental conditions (Harrison and Bates 1991; Groves 
and Lay 1985). Cytotaxonomy seems to be particularly informative, and can be helpful in 
characterizing some species, as the Indian-gazelle, G. bennetti (Furley et al. 1988). Up to 
now, only limited data on the genetic diversity in the genus Gazella are available based on 
protein electrophoresis (Templeton et al. 1987; Granjon et al. 1991). This technique has 
proven to be useful in a number of studies both for captive breeding management purposes 
(see Wayne et al. 1986) and for systematic and phylogenetic studies (see Buth 1984 for a 
review). Advocating a phylogenetic basis for taxonomy, including subspecies groupings 
(see Cracraft 1989 for a discussion of the phylogenetic species concept), the present 
study shows the results of allozyme variation in eight taxa of gazelles and proposes some 
hypotheses on the phylogeny and conservation of this group. 

In particular the following questions were posed: 1. Among the three species present in 
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Fig. 1. Distribution of eight taxa of Gazella in the Arabian peninsula and in East Africa. On Farasan 
Islands lives G. g. farasani. Areas of sympatry for G. g. gazella and G. g. cora 3 G. g. cora and G. g. 
erlangeri are temptative. Data from Groves (1985), Dorst and Dandelot (1972), Harrison and 
Bates (1991) and Kingswood and Kumamoto (1988) 


Saudi Arabia, G. suhgutturosa has been considered very divergent morphologically from 
all other members of the genus, and placed alone in the subgenus Trachelocele (Ellerman 
and Morrison-Scott 1951). Is this morphological divergence accompanied by significant 
molecular genetic divergence? 2. Groves (1983) proposed that the insular gazelles living on 
the Farasan Islands in the Red Sea should be treated as a separate species, namely G. 
arahica. However, Groves (1989) pointed out that G. arahica resembles G. gazella more 
closely than any other species. Based on skull measurements, Thouless and Al Basri 
(1991) proposed to consider gazelles from the Farasan islands as a subspecies of G. gazella , 
named G. g. farasani. Are electrophoretic data able to help differentiate between these two 
hypotheses? 3. At least five subspecies have been attributed to Gazella gazella sensu lato: 
G. g. muscatensis , G. g. gazella , G. g. cora , G. g. farasani and G. g. erlangeri (Groves 
1989; Thouless and Al Basri 1991; Groves et al. 1994) based on morphological 
characteristics. Do some fixed allelic differences characterize any of these taxa, hence 
supporting the concept of separate gene pools in G. gazella sensu lato? 4. The taxonomic 
position of G. mfifrons is not clear. Groves (1975, 1985, 1988) described the Red-fronted 
Gazelle G. rufifrons with 7 subspecies: rufifrons , laevipes, kanuri , tilonura , alhonotata , 
nasalis and thomsoni. The last three subspecies have usually been grouped in a different 
species: G. thomsoni (Dorst and Dandelot 1972; Corbet and Hill 1980; Nowak and 
Paradiso 1983). Is G. thomsoni a subspecies of G. rufifrons ? 

To address these questions, we have analysed by means of protein electrophoresis, 
samples of gazelle species from the Arabian peninsula, and 3 African species. 
































238 


M. Vassart, L. Granjon, A. Greth and F. M. Catzeflis 

Material and methods 

Origin of the animals 

Gazella thomsoni albonotata (n = 8): eight animals from private collections in Saudi Arabia, all 
originating from Sudan. This taxon is called G. rufifrons albonotata by Groves (1985, 1988). 

Gazella dorcas (n = 25): 16 individuals from different private collections in Saudi Arabia, 5 from 
Taif breeding center (National Wildlife Research Center, Saudi Arabia), and 4 from Thumama 
breeding center (King Khaled Wildlife Research Center, Saudi Arabia). The probable origin for all is 
Sudan. These animals would be representative of the taxon G. dorcas isabella mentioned by Groves 
(1969) and Alados (1986/1987). 

Gazella subgutturosa (n = 30): 30 animals from Thumama center. They represent the subspecies G. 
s. marica (Nader 1989). This species is endangered (Iucn 1988) and limited to a very small geographic 
area in the Arabian peninsula, but is still common in Mongolia and in some parts of Russia (Groves 
1988). 

Gazella gazella gazella (n = 16): 16 individuals from the Thumama center. This gazelle is locally 
abundant in Northern Israel. 

Gazella gazella cora (n = 7): 5 animals from Taif center and 2 from a private Saudi collection; they 
may have originated in Saudi Arabia. G. g. cora is rare and endangered in the mountains east of the 
Red Sea (Groves 1988). 

Gazella gazella farasani (n = 5). 4 individuals from a private collection on Farasan Kebir island 
(Red Sea) and 1 wild gazelle from Farasan Kebir island. This taxon, distributed on Farasan islands, 
was previously called G. arabica by Groves (1985) but is more probably a subspecies of G. gazella 
(Thouless and Al Basri 1991). Its population status has been recently examined by Flamand et al. 
(1988), who found that these gazelles were still present in fairly large numbers. 

Gazella gazella erlangeri (n = 15): 6 animals from Taif center, some of them thought to have been 
caught in the South of Saudi Arabia, 5 animals held in private collections in Saudi Arabia, 4 individuals 
from a pet shop in Djeddah. For the latter samples the origin was reported to be the region of Aden 
(Yemen). This subspecies has recently been described by Groves et al. (1994) after morphological and 
skull measurement comparisons with G. g. cora and G. g. muscatensis.. 

Gazella rufifrons (n = 1): 1 specimen from a private collection in Saudi Arabia. The origin of this 
sample is unknown. This gazelle could belong to the subspecies kanuri (synonym = centralis) or 
laevipes (synonym = hasleri ) which are very similar. 


Electrophoresis 

Horizontal starch gel electrophoresis was performed on blood extracts collected by jugular puncture. 
An isotonic saline solution was added to the samples, and the plasma was separated from the red blood 
cells by centrifugation. Samples were duplicated for a reference collection and stored at -30 °C until 
electrophoresis was performed. 

Electrophoresis, staining of the proteins, and scoring the results were conducted according to 
Pasteur et al. (1988). Many allozyme systems were assayed under different electrophoretic condi¬ 
tions, but only 16 loci that gave consistent results were retained (Tab. 1). Statistical and phylogenetic 
treatments were performed without taking hemoglobin into account, as the A and B subunits were not 
separated before running the gels. This prevented us from scoring the different alleles. Nevertheless, 
haemoglobin pattern appears very useful to distinguish G. subgutturosa as will be discussed further. 

For each taxon, the percentage of polymorphic loci (P 95 %, i.e. a locus is considered polymorphic 
when the frequency of the most common allele is not higher than 0.95), mean number of alleles per 
locus (A), and mean heterozygosity (H) were calculated. Genetic distances between samples were 
estimated using formulae from Nei, “unbiased minimum distance” (1978) and Rogers (1972). The 
Arabian oryx ( Oryx leucoryx) was used as an outgroup to root the tree for a phenetic analysis, with 
only 14 loci, because MPI and ACP could not be scored in this species. Allelic frequencies for Arabian 
oryx were published by Vassart et al. (1991). Phenetic analysis was performed via the Distance 
Wagner method (in BIOSYS-1, Swofford and Selander 1989), using the Rogers distance (1972) as 
modified by Wright (1978). 

For cladistic analysis, allelic frequencies were coded as „locus-as-character“ (see Buth 1984). 
Alleles with low frequencies were kept even though this might have introduced a bias for small sample 
sizes. Data for oryx were not used here, as the information they bring after coding was too weak. 
Cladistic analysis was performed using PAUP 3.0 (Swofford 1990). Unordered option was used, and 
the consensus tree obtained with branch-and bound search after 1000 bootstrap replications was 
finally retained (threshold of 95 %). 
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Enzymes 

Tissue 

Locus 

Buffer (pH) 

GOT: Aspartate aminotransferase (EC 2.6.1.1) 

RBC 

1 

TME 6.9/6.9 

ACP: Acid phosphatase (EC 3.1.3.2) 

RBC 

1 

TC 6.4/6.0 

DIA: Diaphorase (EC 1.6.4.3) 

RBC 

1 

TC 6.4/6.0 

ES 10-14: Esterase 10 and 14 (EC 3.1.1.X) 

RBC 

2 

TME 6.9/6.9 

GLO: Glyoxalase (EC 4.4.1.5) 

RBC 

1 

TBE 8.6/8.6 

GPI: Glucose phophate isomerase (EC 5.3.1.9) 

RBC 

1 

TC 6.4/6.0 

LDH: Lactate dehydrogenase (EC 1.1.1.27) 

RBC 

1 

TC 6.4/6.0 

MDH: Malate dehydrogenase (EC 1.1.1.37) 

RBC 

1 

TC 6.4/6.0 

MOD: Malic enzyme (EC 1.1.1.40) 

RBC 

1 

TC 6.4/6.0 

MPI: Mannose phosphate isomerase (EC 5.3.1.8) 

RBC 

1 

TC 6.4/6.0 

NP: Nucleoside phosphorylase (EC 2.4.2.1) 

RBC 

1 

TME 6.9/6.9 

IPO: Superoxide dismutase (EC 1.15.1.1) 

RBC 

1 

TC 6.4/6.0 

ALB: Albumin 

Serum 

1 

LiOH 8.3/8.1 

ES 1: Esterase 1 (EC 3.1.1.X) 

Serum 

1 

LiOH 8.3/8.1 

TRF: Transferrin 

Serum 

1 

LiOH 8.3/8.1 

TME = Tris Maleate, TC = Tris Citrate, TBE = Tris Borate EDTA, LiOH = Lithium Hydroxyde. 


Results 

Among the 8 Gazella taxa, 8 loci were found to be polymorphic; mean heterozygosities 
(H) varied from 0 to 0.085, mean numbers of alleles per locus (A) from 1.0 to 1.31, and 
percentages of polymorphic loci (P) from 0 to 18.7 (Tab. 2). Of the eight polymorphic loci, 
only four exhibited significant variability: ES14, MOD, NP, and TRF (Tab. 2). NP with 
three alleles was polymorphic in five of seven taxa. Rogers (and Nei) genetic distances 
(Tab. 3) varied from 0.012 (0.001) to 0.304 (0.322) among taxa. IPO showed a fixed allelic 
difference between G. thomsoni, G. dorcas and G. rufifrons , on the one hand, and the five 
other taxa, on the other. GOT, DIA and GPI displayed minor variation, with a rare (less 
than 5 % frequency) allele in only one taxon each. 

According to phenetic analysis using the Wagner method (Fig. 2), G. subgutturosa 
clusters together with the 4 samples of G. gazella sensu lato, and appears more similar to 
G. g. erlangeri mainly due to similar allelic frequencies at the NP locus. When hemoglobin 
data are considered, all G. subgutturosa samples had a unique pattern which distinguishes 
that species from all others (Fig. 4). On the phenogram, G. rufifrons lies in an intermediate 
position between G. dorcas/G. thomsoni and the cluster of G. gazella/G. subgutturosa. 

For cladistic analysis, allelic frequencies were coded according to the qualitative method 
described in Buth (1984) (Tab. 4). The unrooted cladogram obtained after parsimony 
analysis (Fig. 3) confirms to a large extent the picture obtained with the phenetic analysis. 
Most of the samples in the G. gazella/G. subgutturosa group are characterized by minor 
autapomorphies (corresponding to low frequency alleles). G. rufifrons and the pair G. 
dorcas/G. thomsoni are clearly distinguishable from each other and from the G. gazella/G. 
subgutturosa group. Two synapomorphies separate G. rufifrons from each of the other two 
clusters. 


Discussion 

These results raise two points that require further discussion: the systematic implications 
of these allozyme data, and their implications for the conservation biology of the Arabian 
gazelles. It should be kept in mind that all the samples used here came from captive 




Table 2. Allelic frequencies for variable loci in 8 samples of Gazella 
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populations, often derived from a few individuals, and 
therefore have probably been subject to sampling ef¬ 
fect and subsequent inbreeding and genetic drift. 
Thus, they may not reflect allele frequencies of the 
parental populations and the results inferred from 
these samples should be taken with caution. Unfortu¬ 
nately, it is difficult to correct for these biases because 
the captive history of these samples is not well 
documented. 

Our allozyme analysis does not support G. subgut- 
turosa as an isolated taxon in the genus Gazella , or the 
concept of the subgenus Trachelocele for this species. 
This is in contradiction to the relationships suggested 
by Ellerman and Morrison-Scott (1951), Hal- 
tenorth (1963), and Groves (1969). Rather, G. sub- 
gutturosa appears to be related to Gazella gazella 
sensu lato. Among the protein systems considered 
here, hemoglobin appears to represent the only dis¬ 
criminative one between these two taxa (see Fig. 4). 

The clustering of G. thomsoni and G. dorcas in the 
cladistic analysis (2 synapomorphies represented by 
ESI4 and MOD) is surprising, since this was never 
suggested by studies of comparative morphology. 
From our data it is difficult to support the view of 
Groves (1975, 1985, 1988) or Lange (1972), who 
described G. thomsoni as a subspecies of G. rufifrons. 
However, chromosome numbers (58) also appear 
identical in G. thomsoni and G. rufifrons (Vassart, 
unpubl. results). Considering allelic frequency (see 
ESI4) G. rufifrons could belong to the G. thomsoni 
group, but with a single specimen for G. rufifrons, this 
conclusion would be speculative. 

There are no standard levels of genetic divergence 
associated with subspecies or species rank (see exam¬ 
ples in Linnell and Cross 1991 and references in 
Cracraft 1989). However, our results suggest that 
several Arabian gazelles do belong to Gazella gazella 
sensu lato, namely, G. g. erlangeri, G. g. farasani, G. 
g. gazella , and G. g. cora. The small genetic distances 
between them are not larger than those found between 
geographic samples in other ungulates (see, for in¬ 
stance, Hartl and Reimoser 1988; Hartl et al. 
1990). The position of G. g. cora well inside the G. 
gazella group does not agree with the hypothesis of G. 
g. cora being the Arabian representative of G. dorcas 
(Groves 1988). This has been confirmed by cytogene¬ 
tic results: G. g. cora, G. g. gazella, G. g. erlangeri and 
G. g. farasani have 35 chromosomes in males and 34 in 
females (V assart et al. 1993) whereas G. dorcas have 
31 chromosomes in males and 30 in females. 

Also, and contrary to the hypothesis of Groves 
(1985) and Groves and Lay (1985), G. arabica (our G. 
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Fig. 2. Dendrogram derived from the Wagner procedure on the matrix of Rogers genetic distances 



Fig. 3. Cladogram obtained with ‘‘locus as character” coding using branch-and-bound search after 
1000 bootstrap replications. Synapomorphies are underlined, other characters are autapomorphies. 
The arrow points to the root leading to the outgroup (see Fig. 2) 

Table 4. Qualitative coding of allelic presence in sample data set from table 2 
See text for explanations 



GOT 

DIA 

ES14 

GPI 

MOD 

NP 

IPO 

TRF 

G. thomsoni 

1 

1 

2 

1 

2 

1 

2 

4 

G. dorcas 

1 

2 

2 

1 

2 

2 

2 

3 

G. rufifrons 

1 

1 

3 

1 

1 

1 

2 

2 

G. subgutturosa 

1 

1 

1 

2 

1 

3 

1 

1 

G. g. gazella 

2 

1 

1 

1 

1 

3 

1 

2 

G. g. cora 

1 

1 

1 

1 

1 

3 

1 

2 

G. g. farasani 

1 

1 

1 

1 

1 

3 

1 

2 

G. g. erlangeri 

1 

1 

1 

1 

1 

4 

1 

2 


g. farasani) should not be treated as a different species, but as a morphological and/or 
geographical race or subspecies of G. gazella: G. g. farasani as proposed by Thouless and 
Al Basri (1991). 

As stated by Furley et al. (1988: page 48), “no common agreement has yet been 
reached on the number of genuine species within [Gazella\. . because of the consider¬ 
able morphological variation shown by some species and the possibility of morphological 
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Fig. 4. Hemoglobin gel with TBE 8.6/8.6 buffer: R; G. subgutturosa , G; G. gazella , T; G. thomsoni . 
G. dorcas has the same patterns as G. gazella 

convergence between taxa as illustrated by the cladistic analysis of skull and jaw characters 
of Alados (1986/1987). For this reason, we feel that a non-morphological approach, such 
as allozyme electrophoresis, should be used for reconstructing the systematics and 
evolutionary relationships of the gazelles. Examination of other species in the genus is 
needed, as well as the scoring of a greater number of loci. 

With respect to conservation issues, Groves (1989) stated that „for conservation 
purposes it is quite clear that each of the Arabian forms [G. subgutturosa marica , G. gazella 
cora and G. dorcas saudiya ] represents a unique gene pool. ,, Our data on the former two 
species support this view. But, as far as local forms are concerned, we have to deal with the 
“dilemma of subspecies” (Ryder 1986) when considering G. g. gazella , G. g. cora , G. g. 
farasani , and G. g. erlangeri. On the basis of our results, they are genetically very similar, 
however, from a conservation biologist’s point of view, they might represent “unique gene 
pools,” each one of which adapted to a particular local environment. 

Among the different guidelines entering into consideration for selecting wild animals to 
be used for breeding purposes, genetic parameters such as mean heterozygosity and 
percentage of polymorphic loci should be taken into account. The data derived from our 
samples of Gazella (see Tab. 2) are similar to those found in natural populations of 
artiodactyls (see Baccus et al. 1983; Vassart et al. 1991, for reviews). The mean 
heterozygosity for 184 species of mammals is 4.1 % (± 3.5 SD) (Nevo et al. 1984), a value 
similar to the one observed for the three Gazella samples represented by more than 10 
individuals in our study. On the other hand, the absence of genetic variability (despite the 
different origins of the samples) for G. g. erlangeri is potentially problematic. This gazelle 
is only known from captive individuals and there is no protected area on its supposed range 
(southwest of the Arabian peninsula). It is possible that the lack of polymorphism and 
heterozygosity resulted from the population bottleneck experienced by this subspecies. 
This could lead to inbreeding problems (O’Brien et al. 1983). 

These results have to be confirmed by other molecular techniques such as mitochondrial 
DNA sequencing. This phylogenetic tool could be useful to differentiate all the different 
species of gazelles even if they are of recent origin, but could be of limited value at the 
subspecific level (Cronin 1992). 
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Zusammenfassung 

Genetische Verwandtschaft einiger Gazella-Arten: Eine Allozym-Untersuchung 

Acht Arten der Gattung Gazella (Bovidae: Artiodactyla) wurden mit einer Allozymanalyse vergli- 
chen, um die genetische Verwandtschaft von Gazellen zu iiberpriifen, die derzeit in Saudi-Arabien zur 
Wiedereinburgerung geziichtet werden. Die elektrophoretische Untersuchung von 16 polymorphen 
Loci deutet darauf hin, daft mehrere ahnliche Taxa, namlich G. gazella gazella , G. gazella erlangeri , 
G. gazella farasani und G. gazella cora zur gleichen Art G. gazella sensu lato gehoren. Vier andere 
Arten sind davon deutlich genetisch verschieden: G. thomsoni und G. dorcas , die genetisch ahnlich 
sind, sowie G. rufifrons und G. subgutturosa. Die Untergattung Trachelocele , der die letzte Art 
aufgrund morphologischer Merkmale zugeordnet worden ist, konnte aufgrund der Allozymvergleiche 
nicht bestatigt werden. Die Polymorphismus- und Heterozygotie-Grade der Gattung Gazella waren 
insgesamt den bei anderen Saugern gefundenen Werten ahnlich. Die Ergebnisse werden im Hinblick 
darauf diskutiert, welche Schluftfolgerungen sich aus den genetischen Untersuchungen fur ein 
angestrebtes Schutzprogramm ergeben. 
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